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Abstract. A version of the Dirac equation is derived from first principles using a 
combination of quaternions and multivariate 4-vectors. The nilpotent form of the 
operators used allows us to derive explicit expressions for the wavefunctions of free 
fermions, vector bosons, scalar bosons, Bose-Einstein condensates, and baryons; 
annihilation, creation and vacuum operators; the quantum field integrals; and C, P and 
T transformations; and to suggest new insights into the meaning of supers ymmetry 
and renormalization. 



1 Introduction 

The algebra required for the quantum mechanical equation for the electron, as 
originally introduced by Dirac, has been discussed by many authors on a formal basis. 
Our aim, here, is to derive it from a more purely physical starting point in a 
fundamental symmetry that we believe exists between the four fundamental 
parameters space, time, mass and charge. This symmetry, we believe, may also be 
significant in explaining many aspects of the Standard Model of particle physics. An 
important feature of the symmetry is that it requires a quaternion representation for 
mass and the three types of charge (electromagnetic, strong and weak), on the same 
basis as we apply a 4-vector representation of space and time, though the three types 
of charge must also be separately conserved, unlike the three dimensions of Euclidean 
space. In applying these two representations simultaneously, we generate a 
quaternion-multivariate-4-vector algebra, which leads to a mathematical 
representation of the Dirac wavefunction as a 4-spinor composed of nilpotent terms. 

In principle, we show that this is the most complete representation of the Dirac 
wavefunction available, as its operators are already quantum field operators, 
eliminating the need for second quantization. We derive explicit versions of the Dirac 
wavefunctions for free fermions, vector bosons, scalar bosons, Bose-Einstein 
condensates, and baryons, together with their intrinsic parities; annihilation, creation 
and vacuum operators; the quantum field integrals; and C, P and T transformations. 
Standard results, such as spin, helicity, the hydrogen atom, and the Schrodinger 
approximation, are easily obtained, while the fact that our operators are also 
supersymmetry operators leads to a completely new insight into the meaning of 
supersymmetry and renormalisation. At the same time, the relation between our 
representation and other standard ones is immediately apparent. These results, we 
believe, arise only because we have started with a formalism which is physically, as 
well as mathematically, true, and leads to a deeper understanding of the physical 
meaning and origin of the Dirac quantum state. 
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2 The Dirac algebra 

Two symmetrical algebras are the basis of the work in this paper. These are the 
algebras of 4- vectors, with real vector units i, j, k and imaginary scalar z; and 
quaternions, with imaginary vector units i, j, k and real scalar 1. The imaginary 
quaternions follow the usual multiplication rules: 

i 2 =f = k 2 = ijk = -\ 
ij = -ji = k 
jk = -kj = i 
ki = -ik = j , 

while the real vector units will be assumed to follow multivariate multiplication rules, 
identical to those for Pauli matrices, and parallel to those for quaternion algebra: 

i 2 = j 2 = k 2 =l 

>j = -j> = ik 
jk = -kj = zi 
ki = -ik = zj . 

In effect, this means defining a 'full product' for two vectors a and b of the form 

ab = a.b + z a x b . 

This extension of the vector system to the full multivariate form will be significant in 
accommodating the concept of spin, and is based on the idea that the two algebras 
used must be completely symmetrical. 

It can be readily shown that the algebra of the gamma matrices used in the Dirac 
equation 

(Y% + im) y/= , 

is completely isomorphic to a combination of the quaternion and 4-vector algebras, 
with assignments of the form 



7° = -ii or 7° = ik 

y l = ik y 1 = ii 

f = jk f=ji 

f = kk f = ki 

f = ij f = ij . 



In general, we will use the first set for the conventional form of the equation, 
obtaining our altered form by multiplying from the left throughout by ij. The cyclic 
nature of the quaternion operators then allows us to redefine our pterins in the form of 
the second set, which will be used for the altered form of the equation. We will, 
however, derive our altered form directly from first principles before relating it to the 
conventional form. 

The complete algebra has 32 parts: 1 real scalar, 1 imaginary scalar, 3 real 
vectors, 3 imaginary vectors, 3 quaternions, 3 imaginary quaternions, 9 real vector 
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quaternions and 9 imaginary vector quaternions. ' ' The existence of 32 parts suggests 
that it can be generated from a binomial combination of five 'primitive' (though 
composite) components, of which the ^matrices are a characteristic set, with terms of 
the opposite sign obtained by reversing the order of multiplication: 

1 

f = ik, f = il f = ij,f = * k , f = V, 

YY = ijh ff = ijh ff = t/X ff = ». ff = -ik, 
YY = ij. YY = **. YY = -»> YY = ikj, ff = ikk, 

^YY = kk, YYY = -kj, y°yY = h Tff = ki,fff = j, 
YYY = k, yW = Y'yY =JX fff = -jj, fff =A 

yoyiyY = j 9 yfff = -/ik, yYyY = »j 5 ffff = -Hi, yYyY = ^. 
Yffff = -i 

Since there is reason to suppose that space, time, mass and charge are the most 
fundamental set of parameters available to physics, 4 ' 5 ' 6 we might expect that the 32- 
part algebra containing their units is of fundamental physical significance, and that 
significant physics will also be contained in the 5-fold structure represented by the 
binomial 'primitives'. 

It is certainly possible, using this algebra, to derive the Dirac equation from the 
relativistic momentum-energy conservation equation 

E 2 — p 1 — m = , 

by, first factorizing and attaching the exponential term e' l(Et ' p ' r) , so that 

(kE + ii p + ij m) (kE + ii p + ij m) e m ' p " r) = , 

and then replacing E and p in the first bracket with the quantum operators, id I dt and 
-i V, to give 

ilc^ + i V+ ijmj (kE + ii p + ij m) e m ' p - r) = . 

This can be written in the form 

( ^ ^ 
ikjr + i V+ ijm y/ = , 

where the wavefunction 

ys= (kE + ii p + ij m) e m " p - r) , 

and the vector elements associated with the second term disappear when we take p in 
a preferred direction, as will often be convenient, (p = l.p can be replaced by p if we 
assume, as is conventional, that only one direction of the vector is well-defined.) With 
the m term fixed as positive, the equation allows for four solutions 
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y/x = (kE + ii p + ij m) e m ' p - r) 
x/h = (kE - ii p + ij m) e m + p - r) 
y/ 3 = (~kE + ii p + ij m) e m + pr) 
y/ 4 = (~kE - ii p + ij m) e m ' pr) 



which represent the four possible combinations of ± E (particle / antiparticle) and + p 
(spin up / down). It can be shown that these solutions are identical to the four 
produced by the conventional Dirac spinor: 

and that the new version of the Dirac equation can be derived from the conventional 
matrix representation and vice versa. 

These are the four solutions as obtained using a common differential operator, 
which may be written: 

( ^ A 
(kE + ii p + ij m) = iJc^ + i V+ ijm 



using the operator interpretation of E and p. However, using a set of four different 
differential operators, 



ilc^ + 1 v+ ijm 



iJc^ + iV+ ijm 



V 



J 



ik^ - 1 V+ ijm 



J 



v 

f d 

ilc^ - i V+ ijm 



V 



which is implicit in the Dirac matrix algebra for the nilpotent wavefunction we have 
defined (see section 4), we may write the equation in matrix form with the solutions: 

y/\ = (kE + ii p + ij m) e i{Et ' p - r) 
y/ 2 = (kE - ii p + ij m) e m " p - r) 
y/ 3 = (-kE + ii p + ij m) e m " p " r) 
y/ 4 = (-kE - ii p + ij m) e' i(Et ' p " r) . 



By using nilpotent wavefunctions, in which the algebra of the differential operator is 
identical to that of the operator part of the wavefunction, we can effectively transfer 
the signs of E and p in the exponential term to the differential matrix, which then 
reduces for each of the individual solutions to the ordinary differential operator. The 
choice becomes a matter of convention, but it is the existence of four differential 
operators, in addition to four solutions, which is responsible for the necessity of 
writing the equation in terms of the product of a 4 x 4 matrix and a 4-component 
vector. 
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The four solutions represent the four combinations of particle and antiparticle, 
and spin up and spin down states. Reversal of the sign of kE produces the 
wavefunction for an antiparticle, while reversal of the sign of zip changes the direction 
of spin. In this form of the equation, the operator p, because it is multivariate and is 
always effectively multiplied by itself in the eigenvalue produced by the differential 
operator, can represent the momentum vector p, the scalar value of p, or the spin term 
o.p. There is no need to multiply the wavefunction operators by any other term, such 
as an additional spinor, to convert momentum states to spin states; spin is built into 
the structure of the multivariate operator p (see, for example, section 13). It is also of 
no significance from a fundamental point of view whether or not we specify a 
preferred direction, as in p x , p y , or p z . We can, in many instances, derive general 
expressions for p from expressions originally derived for such directional terms; the 
equations for the general and the particular are, in each case, identical. As in 
conventional theory, to obtain a scalar probability density ysyf*, we need to use a 
combination of all four solutions. 



3 C-linear maps and lifts 



In our expressions for the Dirac equation and the Dirac wavefunction, p is understood 
to be a multivariate momentum vector with the usual three components. The most 
usual representation for multivariate vectors is the set of Pauli matrices, a , a x , a y , <7 Z , 
which means that we may write the Dirac equation as follows: 



f 3 - 

it^ + ia. V + ijm y/ = , 



v 



J 



where c?. V is understood now to mean the scalar product, o^ x + a y d y + a z d z . This 
has an equivalent 2x2 matrix form: 



kE + ii p z + ij m 



a (p x - ip y ) 



= 0, 



ii (p x + ipy) kE — ii p z + ij m 
which leads to the coupled linear differential equations 

{kE + ii p z + ij m)(j) + ii (p x - ip y )x = , 
ii (p x + ip y )(j) + (kE - ii p z + ij m)x = . 



Choosing the momentum to be along the z-axis (p = p z ), these differential equations 
decouple to leave 

(kE + ii p + ij m)(f) = , 



with 



(kE - ii p + ij m)x = , 

m 



</) = y/\ = (kE + ii p + ij m) 



-KEt-pz) 
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X= ¥2 = (kE - ii p + ij m) 



-i(Et + pz.) 



These are the positive energy solutions; the quaternion representation allows a simple 
deduction of the negative energy solutions: 



y/3 = (-kE + ii p + ij m) 



\j/ A = (-kE - ii p + ij m) 



i(Et + pz) 



i(Et - pz) 



To carry out quantum mechanical calculations, we need to define a scalar product for 
the wavefunctions. However, the wavefunctions in the quaternion / spinor form are 
effectively a tensor product of two representations. 

It is convenient to use a C-linear map F: tf — > C 2 is a mapping from complex 
quaternions, tf (q^ M , q M e C) to C 2 . 



F: if -> C 2 , Q -> Q 
where Q is a complex quaternion, and 
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Applying this to ^i, we obtain 



f _/ ' 








® 








v0. 



where K = p I E, £ = m I E, and ® is the tensor product between the two 
representations. Hence all four solutions may take the form 



V2 



f ~i \ 


K+ i£ 
\ J 



f \ 

— i 


\ v -/r+ z£ y 



^3 = 





K+ i£ 

\ J 



W4 = 



V"2 



z' \ 

i 





This mapping leads to the conventional Dirac equation for these states, which in the 
case of y/\ is 
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f -iE 


im + p 





-iE 


im - p 



-im + p 

/£ 




\ 

-im - p 


iE J 



( -i \ 


k+ is 



= o 



Taking the four solutions, y/\, ifr, ifi, y^, as defined above, we combine their values 
of 

y/ y/= y/* \j/ , 

to obtain the probability density 

4(-E 2 -p 2 -m 2 ) = -SE 2 , 

which becomes 1 on application of the normalising factor i I sJlE for each solution. 
y/*, of course, defines a row vector with the terms, y/\*, ifi*, ife*, ifa*, while ^is the 
column vector with y/\, ife W3, Ya- 

For a Dirac wavefunction of the form 

y/= (kE + ii p + ij m) e i{Et ' p " r) 

we have also defined an Hermitian conjugate, 

yj = fy/f =yA = (kE- ii p - ij m) e m ' p ' r) 

and an adjoint wavefunction, 

y/= y/f = y/ ik = (- iE -7'p + i m) e i(Et ' p r) , 

though it may be equally convenient to replace this with the same function multiplied 
from the left by ik, so that 

yr=(kE + iip + ij m) e i(E, - p - r) , 

as this function equally satisfies the adjoint equation 

( - ^ 



-^■ik + Vyd- ynjm 



= 



and is only ever used as a multiplier from the left. 

From this, we may derive the bilinear covariants, and hence the current density 
terms (using the four solution sum and the normalising factor): 

yfy/= ynky/= ik E 2 IE 2 = ik 
yr/ [ y/= yAi yr= -ii p 2 IE 2 
yrfyr= yAjy/=-ijp 2 /E 2 
yfy/= yAky/=-ikp 2 /E 2 
yrfy/= ytijy/= - ij m 2 IE 2 . 
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The first four quantities (the scalar and vector terms) are the components of the 
current-probability density 4-vector, F = y/fys, where cPF = in the absence of 
external fields. Using the original form of the adjoint wavefunction, the scalar, vector 
and pseudoscalar terms become: 

ysfy/= 1 

Wfy/=-ijjp 2 /E 2 
yrfy/=-ijkp 2 /E 2 
yrfy/=i mIE 2 , 

while the pseudovector terms take the values: 

ysffys=-jp 2 /E 2 
i/rfy 1 y/=-i 

yrff ys=-k, 

and the tensor terms are zeroed. 

We may also use the bilinear covariants to construct the zero Dirac Lagrangian, 
which results when the equations of motion are obeyed: 

L = y/^d^ - m) y/= . 

4 The quaternion form derived from a matrix representation 

We have shown, in the previous section, that our quaternion nilpotent algebra can be 
used to generate the standard version of the Dirac spinors. The reverse process is also 
possible, but requires a more subtle argument which suggests why exactly four 
solutions are required (only four are, of course, possible if they are to be 
independent). It also suggests a link may be found with Dirac algebras that structure 
the four solutions as components of a four- vector or quaternion set. 

Here we interpret the matrices as dyadics formed from quaternion (or, 
alternatively, 4-vector) components, arranged by row and column. Starting with the 
equation 

(a.p + /3m -E) y/ = 
and taking (for convenience, without loss of generality) p = p y , we obtain 





( 
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Oty = 
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Using 
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and applying the unit 4x4 matrix to E, the Dirac equation becomes 




f -E 





im 




-ip \ 









-E 


ip 




im 


^2 


(oc.p + f5m -E) y/ = 


—im 


-ip 


-E 







^3 




V ip 


-im 
















where the column vector is the usual 4-component spinor, and the terms E and p are 
the quantum operators which give the eigenvalues represented by these symbols when 
applied to the exponential term of the wavefunction. (It will be convenient here to 
refer to the components of the 4x4 matrix in terms of these eigenvalues rather than 
in terms of the operators which produce them.) 

We can interpret the rows and columns of this matrix as having either 4-vector or 
quaternion coefficients. Let us choose the quaternion operators 

j, i, k, 1 

as the respective coefficients of the 4 rows. The 4x4 matrix now becomes a single 
row bra matrix with the columns: 

-jE - ikm + ip 

— iE — ikp — im 
ijm + Up — kE 

- ijp + Urn - E . 

If we multiply these terms from the left by the respective column coefficients 

i, ~j, l,k, 

we obtain in each case the expression 

-kE + ii p + ij m , 

which, when multiplied from the right by a wavefunction beginning with this term, 
gives a zero product. In order to show that this is equivalent to the Dirac equation in 
matrix form, we need to show that the four solutions y/\_, yti, ifi, y/<\, multiplied by the 
appropriate quaternion row coefficients, each result in expressions beginning with this 
term. 
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Suppose, therefore, that we have the four solutions, as previously assumed: 



y/ x = (kE + ii p + ij m) e iiEt - p - r) 
y/ 2 = (kE - ii p + ij m) e m " pr) 
y/ 3 = (-kE + ii p + ij m) e m " pr) 
y/ 4 = (-kE - ii p + ij m) e i{Et - p - r) 



We can show that the terms 



k y/ x k 

-j WJ 

1 y/ 3 1 

i y/^i 



each produce the expression 



(-kE + iiv + ijm)e- i{E '- p - r) . 
So, multiplying each of the terms 

k y/ x 

-j ¥2 
1 ¥ 
i y/ 4 

from the left by (-kE + ii p + ij m) results in a zero product. 

For convenience, we may rearrange these to give a ket matrix of the form: 



y/2 

W3 



= i y/ 4 -j y/ 2 + l yA + ky/y 



where the row coefficients are identical to the column coefficients of the bra matrix 
(and even allow the elimination of the - sign before -j). The resulting equation is 
identical to the quaternionic Dirac equation, which we have previously derived by 
direct means, with the four solutions representing the four possible combinations of + 
E and + p states. 

The derivation demonstrates that the reason for the use of 4 x 4 matrices is, in 
fact, the fundamentally quaternionic nature of the Dirac wavefunction. Ultimately, 
this is because of the 4-vector space-time used in relativistic equations, which is 
symmetrical with the quaternion algebra for mass and charge. In our understanding, 
the use of quaternionic operators to define the weak, strong and electromagnetic 
charges (w, s, e) maps directly on to the use of the same operators for the terms E, p, 
m in the wavefunction, the existence of these terms as independent units stemming 
ultimately from the separate identities of the fundamental parameters time, space and 
mass (T, S, M). Because quaternion operators define the meaning of the rows and 
columns in the Dirac matrix, the only way we can map charges w, s, e on to the E, p, 
m terms via these operators is to have 4x4 matrices, and hence a 4-dimensional 
space-time signature in the equation. 
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It is significant that the application of quaternion operators to w, s, e and E, p, m 
- and, by implication, to the more fundamental parameters time, space and mass (or 
energy) (T, S, M) - is incomplete, each of these groups of three terms requiring a 
fourth to complete it. We can identify the respective fourth terms, without difficulty, 
as mass (M), angular momentum (J) and charge (Q). The 4x4 Dirac matrix, in effect, 
incorporates the fourth term as a zero quantity. But, while the matrix requires four 
quantities to which the columns and rows apply, the 4-component spinor allows only 
four possible solutions from their combination. Interpreting the solutions in terms of 
the number of relative sign combinations of the component terms allows only three of 
the terms to be nonzero. Significantly, the excluded term in each case is an invariant, 
the system requiring only one invariant quantity (e.g. m or M) to demonstrate the 
variability of the others (E, p, or T, S). 

5 Some aspects of the Dirac algebra 

The power of the current version of the Dirac algebra lies in the fact that it uses 
nilpotents (or square roots of zero) rather than more conventional forms of the Dirac 
wavefunction. The two may, however, be easily related. Consider the equation: 

(iE - ik p - im) (iE -ik\) + im) = , (1) 

which is just one version of the usual energy-mass-momentum equation. The two 
bracketed terms are clearly different, and so are not square roots of zero or nilpotents. 
However, if we multiply from the left by -j and from the right by j (in the process 
changing the Clifford algebra of the operators from CLy to CI23), we get: 

-j(iE - ik p - im) (iE - ik p + im) 7 = 0. 

This now becomes: 

(kE + ii p + ij m) (kE + ii p + ij m) = . (2) 

The two bracketed terms are now identical, and so each is a nilpotent or square root of 
zero. Equation (1) is, in effect, the usual way of writing the Dirac equation, with the 
left-hand bracket becoming the differential operator and the right-hand bracket the 
wavefunction, whereas (2) is just as valid, but much more powerful. In incorporating 
the m term directly, it effectively goes beyond the conventional Dirac wavefunction 
towards a quantum field interpretation, though it retains, as we have shown, all the 
physical interpretation available to the conventional term. 

It is important to stress here that what is presented as a purely formal 
development has a physical, rather than mathematical, origin. It is based on the idea 
that the four fundamental parameters of physics, space, time, mass and charge, have 
fully symmetrical properties, one aspect of this being that, while space and time can 
be represented by a 4-vector, mass and charge have a symmetrical formal 
representation as a quaternion, the real part of which is represented by mass, and the 
three imaginary parts by the sources of the electromagnetic, strong and weak 
interactions (e, s, w). 4 ' 5 ' 6 The development of multivariate approaches to vector 
algebra by Hestenes et al. 7 suggests, in addition, that the vector part of the 4-vector 
space-time is multivariate or quaternion-like, manifesting its additional properties 
through spin. It then seems reasonable to suppose that what appears to be the most 
fundamental equation for describing particle states should be described by an algebra 
combining fully symmetrical 4- vectors and quaternions. 
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Though there are similarities in some aspects to the quaternion and double 
quaternion algebras used by some authors, 9 " 17 the differences are significant in 
creating a version of the Dirac formalism based on nilpotent operators, in particular 
the extra element of complexification which this requires. The use of vectors, in 
addition to quaternions, certainly helps to distinguish the two algebraic systems 
needed, but they are used here because they are believed to be a true physical 
representation of the components involved, and, at the end of the paper, they will be 
required for a physical explanation of the Dirac state. The complexification is 
physically relevant because it comes from the complexification of the E term (or, 
equivalently, of time), as it does in classical relativistic physics. 

The physical connotations of the quaternions and vectors also explain the choice 
of symbols used, though some readers may find it occasionally difficult to distinguish 
the vectors i, j, k from the quaternions i, j, k. In fact, for the most part, the vector 
aspect is subsumed in the p term, and we only need the quaternion symbols, retention 
of the vectors being optional. In addition, vectors and quaternions can be regarded, in 
some sense, as 'numbers', rather than as merely complicated abstract mathematical 
objects. This means that we can write down precise expressions for the wavefunctions 
which have an immediate physical meaning - a meaning which comes from the fact 
that they originate in a mathematical representation of space-time and mass-charge 
which is itself physically meaningful. 

It is significant that, in using the mathematical formalism which results from this 
physical origin, it becomes apparent that what the Dirac equation set out to do in 1928 
can be achieved more conveniently by using variations from the particular form used 
at that time. In principle, the Dirac equation, as used by Dirac, requires wavefunctions 
to be ideals, incorporating only a restricted part of the algebra. The above analysis 
shows, however, that limiting the wavefunction to an ideal limits drastically the 
physical interpretation that one can derive from the algebra. In other words, a 
restriction of the wavefunction to part of the algebra is a consequence of the fact that 
the traditional form of the Dirac equation is incomplete - a fact which results in the 
necessity of second quantization. We will show, however, in this paper, that second 
quantization is unnecessary. The Dirac operators as used here are already quantum 
field operators, and it is the multiplication of the wavefunction by j or i'f to create the 
nilpotent form that makes this possible. It can also be shown that the operators are, in 
addition, annihilation and creation operators, vacuum operators, and supersymmetry 
operators, and that they can explain fermion and boson wavefunctions, spin Vi, Pauli 
exclusion, CPT, etc, and even, it would seem, renormalization. 

Now, the 'standard' form of the Dirac wavefunction, as we have defined it, is 
eigenpotent, when the wavefunction is written in the form it has before 
postmultiplication by j, and each of the four terms in the column vector is an 
idempotent. One of the strengths of the nilpotent version here defined is that, both the 
column vector and each of its individual terms is a nilpotent, which is the extreme 
version of both eigenpotent and idempotent. A minimal left ideal is a column vector 
of 'primitive' idempotents multiplied from the left, which is a subgroup of some 
larger group. The (kE + zip + ijm) form is a column vector, multiplied from the left, 
composed of idempotents before multiplication by j, and nilpotents (an extreme form 
of idempotent) after multiplication, and it is, in effect, a subgroup of a larger group. It 
is an extreme case of an ideal. 
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The five objects represented by the gammas do not include 1 as such because 
there are no strictly scalar terms in the Dirac equation in its nilpotent form. The more 
standard form of Dirac has i (complex number) or 1 for the mass term but no f. The 
nilpotent formulation effectively replaces 1 or i with f . In the pentads from the 
double quaternion or vector-quaternion set (say i, j, k, i, j, k), each pentad must 

1 g 

incorporate one (and only one) triad such as ii, ij, ik, or equivalent. This is important 
physically because it 'privileges' the vector nature of p. The nilpotent Dirac algebra 
works for p as a vector in 3-D, for p as a scalar value, or for a.p in a defined direction 
as required, and the formulation also works for E and p incorporating field terms, or 
producing covariant derivative equivalents (with, of course, the corresponding 
changes made in the values of E and p, e.g. E = dldt - ie(j) and p = -iV + eA). 

6 Fermion and Boson Wavefunctions 

The vector-quaternion Dirac equation is simpler and more powerful than the 
conventional representation, and allows a more direct insight into the equation's 
physical meaning. In the first place, we have a direct expression involving E, p and m 
for the wavefunction, and, secondly, we can immediately see the explanation for such 
things as fermionic- and bosonic-type wavefunctions. Since the exponential terms in 
the wavefunctions multiply as scalars, and can be expressed, in the nilpotent algebra, 
in a common form for all wavefunctions, we will often find it convenient, here, to 
discuss the quaternion operators such as (kE + ii p + ij m) as the 'wavefunctions' 
when we are examining the properties of superposed states. 

The immediately obvious aspect of fermionic wavefunctions is that they are 
composed of noncommutative operators, and hence are clearly antisymmetric. Also, 
multiplying two identical fermion wavefunctions produces zero, as in 

(kE + ii p + ij m) (kE + ii p + ij m) = . 

As nilpotents, or square roots of zero, the wavefunctions have many advantages over 
conventional representations. Pauli exclusion is immediately obvious because the 
product of two identical wavefunctions y/ys, each with four solutions, is zero because 
each term in it is zero. That is, 

(kE + ii p + ij m) (kE + ii p + ij m) = ; 
(kE - ii p + ij m) (kE - ii p + ij m) = ; 
(-kE + ii p + ij m) (-kE + ii p + ij m) = ; 
(-kE - ii p + ij m) (-kE - ii p + ij m) = . 

Reversing the sign of either kE or ii p or both in one of the factors, however, 

2 2 2 2 

always produces a nonzero scalar product, such as a multiple of E +p - m , or —E + 

2 2 

p — m , when summed up over the four solutions. Bosons are superposed states of 
fermion and antifermion, so are a product of terms with opposite signs of E in the 
components. A vector boson (spin = 1) has fermion and antifermion with the same 
sign of p. So, its four component terms are: 

(kE + ii p + ij m) (-kE + ii p + ij m) ; 
(kE - ii p + ij m) (-kE - ii p + ij m) ; 
(-kE + ii p + ij m) (kE + ii p + ij m) ; 
(-kE - ii p + ij m) (kE - ii p + ij m) . 
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2 2 2 2 

The sum is a nonzero scalar, 4 (E + p + m ) = 8 E , before normalization. This is 
also true if the meson is massless, because the sum of the terms: 

(kE + a p) (-kE + ii p) = ; 
(JcE - ii p) (-kE - ii p) = ; 
(-kE + ii p) (kE + ii p) = ; 
(-kE - ii p) (kE - ii p) = 

2 2 2 

is still nonzero, 4 (E + p ) = 8 E , before normalization. This means that vector 
bosons can be massless (e.g. gluons, photons). For a scalar boson (spin = 0) the p 
terms of fermion and antifermion have opposite signs. Hence the product becomes the 
sum of: 

(kE + ii p + ij m) (-kE - ii p + ij m) ; 
(kE - ii p + ij m) (-kE + ii p + ij m) ; 
(-kE + ii p + ij m) (kE - ii p + ij m) ; 
(-kE - ii p + ij m) (kE + ii p + ij m) . 

2 2 2 2 

This is, again, a scalar, 4 (E - p + m ) = 8 m , before normalization. However, if we 
now make the mass zero, we then obtain the sum of: 

(kE + ii p) (-kE - ii p) ; 
(kE - ii p) (-kE + ii p) ; 
(-kE + ii p) (kE - ii p) ; 
(-kE - ii p) (kE + a p) . 

2 2 

This is zero for a massless particle, 4 (E - p ), and, hence, no scalar boson can have 
zero mass. If the Higgs, for example, is a scalar boson, then it must have mass. So 
must the pion, while the massless gauge particles (gluons and photons) must all have 
spin 1. A boson with a scalar wavefunction will, in addition, be incapable of 
description via the Dirac equation, as the eigenvalue of a quaternionic quantum 
operator can only produce a zero product when multiplied by another quaternionic 
term. 

There is one other possibility. If we produce a composite wavefunction, of which 
both parts are fermions (or antifermions), but spins are opposite, then we get a product 
which is the sum of: 

(kE + ii p + ij m) (kE - ii p + ij m) ; 
(kE - ii p + ij m) (kE + ii p + ij m) ; 
(-kE + ii p + ij m) (-kE - ii p + ij m) ; 
(-kE - ii p + ij m) (-kE + ii p + ij m) . 

2 2 2 2 

This is a nonzero scalar, 4 (-E - p + m) = -8 p , before normalization, and so 
theoretically can exist, and this is what happens with a Bose-Einstein condensation. 

Baryon wavefunctions may also be constructed from the nilpotents. The standard 
QCD representation is the antisymmetric colour singlet of SU(3): 

y/~ (BGR - BRG + GRB - GBR + RBG - RGB). 
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Here, we use a mapping such as: 



BGR (kE + ij m) (kE + ij m) (kE + ii p+ ij m) 

- BRG (kE + ij m) (kE - ii p + ij m) (kE + ij m) 
GRB (kE + ij m) (kE + ii p + ij m) (kE + ij m) 

- GBR (kE + ij m) (kE + ij m) (kE - ii p + ij m) 
RBG (kE + ii p + ij m) (kE + ij m) (kE + ij m) 

- RGB (kE - ii p + ij m) (kE + ij m) (kE + ij m) , 

2 2 

with each term equivalent to -p (kE + ii p + ij m) or -p (kE - ii p + ij m). 

It is noticeable here that there are three cyclic and three anticyclic combinations, 
just as in the conventional representation. Elsewhere, 19 ' 20 we have represented the 
colour exchange process, by the exchange of the p term, as in conventional gluon 
theory. This fits in with p and strong charge being isomorphic, both being three- 
dimensional, but with only one direction well-defined. Because there is only one spin 
term, it also predicts that the spin is a property of the baryon wavefunction as a whole, 
not of component quark wavefunctions. This structure is determined solely by the 
nilpotent nature of the fermion wavefunction. Put in an extra p into the brackets 
missing them, and we immediately reduce to zero. 

With the spinor terms included, each of these is represented by a tensor product 
of three spinors, for example: 



(kE + ij m) (kE + ij m) (kE + ii p+ ij m) 



where 



® 



® 



e 




So this representation encompasses both spin Vi and spin 3/2 baryon states. 
7 CPT Symmetry 

The quaternionic version of the wavefunction also clarifies the meaning of CPT 
symmetry. Each of the quaternion operators is responsible for one of the three 
component symmetries. C, P and T operations on a quaternion wavefunction may be 
represented as the following operations: 



P 
T 

C 



i (kE + ii p + ij m) i = (kE - ii p + ij m) 
k (kE + ii p + ij m) k = (- kE + ii p + ij m) 
-j (kE + ii p + ij m) j = (-kE - ii p + ij m) . 



From this we may see that: 



CP = T 
PT = C 
TC = P 



-j (i (kE + ii p + ij m) i) j = k (kE + ii p + ij m) k = (-kE + ii p + ij m) 
i (k (kE + ii p + ij m) k) i = -j (kE + ii p + ij m) j = (-kE - ii p + ij m) 
k (-j (kE + ii p + ij m) j) k = i (kE + ii p + ij m) i = (kE - ii p + ij m) 



and that TCP = identity, since: 

k (-j (i (kE + ii p + ij m) i) j) k = -kji (kE + ii p + ij m) ijk = (kE + ii p + ij m) 
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A violation of charge conjugation, such as happens (at least partially) in the weak 
interaction, would effectively make the ij m term unable to assume the opposite sign, 
and so would lead to something like 



j YJ =(k E + ii p + ij m)= yr. 

Similar arguments could be applied to parity and time -reversal violation. In such a 
case, since 

j YJ = (k E + ii p - ij m) = yr- 2ij m , 

that part of the wavefunction involving the weak interaction would require a term like 
2ij m to be equivalent to 0. For a neutrino, which only has a weak component of 
charge, and involves violation of charge conjugation for the entire wavefunction, this 
would require, in the first instance, a particle with zero mass. 

In general, to reduce a term like + 2ij m to requires adding its negative value 
elsewhere. Weak interactions require such additions for fermions where violation of 
charge conjugation makes the weak interaction unable to distinguish the weak charge 
components +w and -w. 

8 Parities of bosons and baryons 

Defining the parity transformation on ^as i y/i, we can now investigate the intrinsic 
parities of ground state bosons and baryons. Applying the transformation to a scalar 
boson, we obtain: 

i (kE + ii p + ij m) {-kE - ii p + ij m) i ; 
i (kE - ii p + ij m) (-kE + ii p + ij m) i ; 
i (-kE + ii p + ij m) (kE - ii p + ij m) i; 
i (-kE - ii p + ij m) (kE + ii p + ij m) i. 

If we take the first term, and use -ii= 1, then the parity transformation produces: 

—i (kE + ii p + ij m) i i (-kE - ii p + ij m) i . 

Now, we have a parity transformation on each bracket, with an additional - sign. This 
produces: 

- (kE - ii p + ij m) (-kE + ii p + ij m) . 
Applying to each of the terms, we obtain: 

- (kE - ii p + ij m) (-kE + ii p + ij m) ; 

- (kE + ii p + ij m) (-kE - ii p + ij m) ; 

- (-kE - ii p + ij m) (kE + ii p + ij m) ; 

- (-kE + ii p + ij m) (kE - ii p + ij m) . 

The total transformed wavefunction iy/i thus becomes -yr. The original wavefunction 
therefore has negative parity. 
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For the vector meson: 

- (kE - ii p + ij m) (-kE - ii p + ij m) ; 

- (kE + ii p + ij m) (-kE + ii p + ij m) ; 

- (-kE - ii p + ij m) (kE - ii p + ij m) ; 

- (-kE + ii p + ij m) (kE + ii p + ij m) . 

Again, i y/i becomes - y/, and the original wavefunction has negative parity. 
Let us try the same operation on a baryon. Taking one of the terms: 

(kE + ij m) (kE + ij m) (kE + ii p+ ij m) , 
we apply a parity transformation: 

i (kE + ij m) (kE + ij m) (kE + ii p+ ij m) i . 
This time, we can write it in the form: 

i (kE + ij m) i i (kE + ij m) i i (kE + ii p+ ij m) i , 
with no sign change. This term becomes: 

(kE + ij m) (kE + ij m) (kE - ii p+ ij m) . 

Taken over all the terms (three with p, and three with - p), then: 

(kE + ij m) (kE + ij m) (kE - ii p+ ij m) 
(kE + ij m) (kE + ii p + ij m) (kE + ij m) 
(kE + ij m) (kE - ii p + ij m) (kE + ij m) 
(kE + ij m) (kE + ij m) (kE + ii p + ij m) 
(kE - ii p + ij m) (kE + ij m) (kE + ij m) 
(kE + ii p + ij m) (kE + ij m) (kE + ij m) , 

and 

i y/i = y/. 

The baryon wavefunction has positive parity. 

These calculations, of course, apply to the ground state values only, because if 
extra angular momentum terms are added, then extra terms must be supplied to the 
wavefunctions, the sign of parity reversing for each additional term. 

9 Annihilation and creation operators 

Yet another significant aspect of the quaternion Dirac algebra using nilpotents is that 
it is already effectively second quantized. We have seen that the quaternion operator 
(ikd Idt + i V+ ij m) acting on a state 

(kE + iiv + ijm) e - i(Et - p - r) 

leads to the equation: 

(ikd Idt + i V+ ij m) (kE + ii p + ij m) e~ m ~ p - r) = 0. 
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This operation with operator S = (ikd/d t + iV+ ij m) may be thought of as a creation 
operation acting on the single particle fermion state which is already filled. The result 
is therefore zero. We may obtain the corresponding annihilation operation by finding 
the Hermitian conjugate of (ikd/d t + i V+ ij m). 
We know that: 

S e - iiEt - p - r) = (kE + ii p + ij m) e - i{E '- p - r) . (1) 

3 is acting as though it were the creation operator at acting on the group of 
translations and rotations that we call vacuum. We can therefore write at as 

af = (V2E) (kE + ii p + ij m) . 

Now, the Hermitian conjugate of expression (1) is: 

P e -m - P .r) = ( _ kE + iiv + ijm) e m - P .r) } (2) 

Here, because it describes the creation of the anti-particle, J f is acting as though it 
were the annihilation operator a conjugate to at, so that 

a = (1/2E) (-kE + ii p + ij m) . 

The operators have been modified by a factor (l/2£) so that, when acting on a state, 
ata will reproduce the state up to a scalar multiplicative factor which does not affect 
the quaternion properties of the state. 

It is easy to verify that these two operators have the commutation relations 
appropriate to fermion annihilation and creation operators: 

aat + ata = 1 . 

We now need to find the vacuum state that a 1 can act upon to lead to the single 
particle state: 

(kE + iiv + ijm) e - i{Et - p - r) . 

We need only consider the following: 

a (kE + ii p + ij m) e~ i{Et ~ p - r) = VAC 
= IE (V2E) (E -ijp + ii m) e^''^ . 

To check this we need only consider 

at (a (kE + ii p + ij m) e~ m ~ p - r) = a + VAC 
= at (V2E) 2E (E -ijp + ii m) e m ~ pr) 
= (1/2E) 2E (kE + ii p + ij m) e~ m ~ p - r) 
= (kE+iip + ijm) e - i{Et - p - r) . 

It is easy to verify that the further action of a onto VAC leads to zero. 
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10 The quantum field 

Let us define the following: 

The creation operator for an electron, spin up: 

a|(p) = {kE + zip + ijm) I IE 
The annihilation operator for an electron, spin up: 

fl(P) = {-kE + zip + ijm) I 2E 
The annihilation operator for an electron, spin down: 

at(-p) = (kE - z'ip + ijm) I IE 
The creation operator for electron, spin up: 

a(-p) = {-kE - z'ip + ijm) I 2E 
The annihilation operator for a positron, spin down: 

b(p) = {kE + z'ip + ijm) 1 2E 
The creation operator for a positron, spin down: 

&t(p) = {-kE + z'ip + ijm) I IE 
The annihilation operator for a positron, spin up: 

b{-p) = {kE - z'ip + ijm) I IE 
The creation operator for a positron, spin up: 

b\ (-p) = {-kE - z'ip + ijm) I IE 
Now, the anticommutators, 

{af(pi),a(p 2 )}= ^(pi-p 2 ) 

and 

{fct(Pi),a(P2)}= <?(Pi-P2) 

are valid for 

(£ 1 -£ 2 ) 2 -(p 1 +p 2 ) 2 -(m + m) 2 = . 



Writing out the Fourier superpositions of all possible states, we have 











j(fl(p) 




+ fl(-p) 





+ (fct(p) 











1°; 


+fct(-p) 


v 1 ; 







x) = JVp{(at(p)(l 0) + at(-p)(0 l))c*"+(*(p)(l 0) + K"P)(0 1))<T**} 



Here, because of the explicit expressions used for the creation and annihilation 
operators, we need only use 1 and states in the spinors. 
Now 

yKx) = I E\ p> , 

where 

af (p) | 0> = I E, p> 
at(-p) I 0> = \E, -p> 

fct(p)l0>= |-£,p> 
fct(-p)l0>= |-E,-p> 



^(x)|F,p') = Jj 3 p (a(p) 













l(fl(p) 




+ a(-p) 









^1 \ 



+ *>t(-p) 



V ! 7 



)e m \E\\ 



Since 



therefore 



a(p) | p'> = a(p) af(p) I 0> 

= (<?( Pl -p 2 )-a(p) af(p)|0> 
= |0> , 



(0\y(x)\E\p) = 



In principle, therefore, we no longer need to define an explicit process of second 
quantization. Our nilpotent operators are already second quantized. 

11 Spin 

In addition to such radical reformulations, various standard results can be 
accommodated into the new formalism by replacing the gamma matrices with 
quaternion operators. Particularly important is the derivation of fermion spin. In the 
conventional treatment of spin, we write 
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Also, 

and 

while 



So, 



or 



and 



[o, = [a, iyj.p + %m] . 
&i = iW57i . with/= 1,2,3 



ri = m 

^ = ki 

75 = • 



A 

£7i = -1 

A 

02 = -J 

<7 3 = -k 



o = -l , 



y = il , 



where 1 is the unit (spin) vector. 

Since fom = ikm has no vector term and o no quaternion, they commute, and we 
may derive the conventional 

[a, Yom] = 

and 



Now, 



So, 



[a, X] = [a, i ^y.p] 



('WP = -7 (ipi + + kp 3 ) 



[a, = 2y (ijp 2 + ikp 3 + J>i+ jkp 3 + kipi + kjp 2 ) 
= 2ij (k(p2 - Pi) + j(pi -P3) + i(p3 -P2» 

= 2// 1 x p . (3) 

In more conventional terms, 

[a, #] = 2*i (k0? 2 - Pi) + j(pi - pz) + i(p 3 - p 2 )) 
= 2/A: y x p 

= 2x,yxp. (4) 

Simultaneously, if L is the orbital angular momentum r x p, 

[L, H\ = [r x p, z'}{,Y.p + x,m] 
= [r x p, z'^Y-p] . 
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Taking out common factors, 



[L, ti\ = i # [r, y.p] x p 
= -ki [r, l.p] x p 
= -j [r, l.p] x p 



Now, 



[r, l.p] y/= - ii 



dy/ djxy/) 
dx dx 

= il y/. 



y 



dy/_ J(yy) 
dy dy 



ik 



dy/ djzw) 
dz dz 



Hence, 



[L, «]=-jflxp. 
This, again, can be converted into conventional terms: 
[L, ti\ = iki 1 x p 

= i»yxp- 

Using either (3) and (5) or (4) and (6), we may write 

[L-l/2,#] = 
[L + 8/ 2, = . 



(5) 



(6) 



or 



Hence, (L - 1 / 2) or (L + a/ 2) is a constant of the motion. The physical origin of 
spin states will be discussed below. 

12 Helicity 

The term 

o.p = -pi -pi-P3 = -p 

is defined as helicity, and, since it has no vector or quaternion terms, and has only 
terms of the form dldx, dldy, and dldz in common with 

i'M-P = -J + SP2 + kp 3 ) 
and also clearly commutes with y^m = ikm, then 

[o.p,#]=0 

and the helicity is a constant of the motion. 

For a particle with zero mass, the term k E + ii p + zy m reduces to k E + a />, 
where p actually represents o.p. E also becomes equal to + p. For positive energy 
states, 

£ = o.p . 
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So the spin is aligned antiparallel to the momentum (has left-handed helicity). Then, 

ij (k E + ii p) = ij (k - ii ) E = (ii -k) E 
and the spinor wavefunction follows the rule: 

ij u L = -u L . (7) 

For negative energy states, 

E = - a.p , 

In this case, the spin is aligned parallel to the momentum (has right-handed helicity). 
Then, 

ij (k E + ii p) = ij (k + ii ) E = (ii + k) E 
and the spinor wavefunction follows the rule: 

ij u R = u R . (8) 



From (7) and (8), we may derive the relations 



and 



r i - a 

{ 2 j 



2 



Ur = 



V 2 J 



Ur — Ul 



Ur = 



If we define the right- and left-handed components of the wavefunction, y/ R and y/ L by 
the expressions 



and 



1 + 75 

V 2 J 



¥= 



¥-- 



(I + ij 



a - v 

2 



¥, 



we may derive the explicit expressions 



and 



f \ + iL ] 
2 



(k E + ii p + ij m) = 



1 - ij 

— Uk E + ii p + ij m) = 



1 + H 

v 2 J 



(k (E + p) + m) 



-km 



j^\(k(E + p) + m) 



(9) 
(10) 



It is apparent from (9) that y/R = when m = and E = -p = a .p. 
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13 The Schrodinger approximation 

The derivation of the Schrodinger approximation is purely conventional, and doesn't 
require the specific use of the nilpotent wavefunction. We take the Dirac equation in 
the form 

O'y.p + m - pE) y/ = 

and choose, again without loss of generality, the momentum direction ip x = p. Here 
again, also, E and p represent the quantum differential operators, rather than their 
eigenvalues. This time, we make the conventional choices for /?: 







r\ 










^ 









1 






















-1 





and for 




V° 




























f° 
























-i 












i 














V 










°y 


leading to the representation: 












/ E — m 












-p 




E 


- m 




-P 












P 




E- 


m 







I P 












-E - m 



¥3 



= 



This can be reduced to the coupled equations: 

(E-m)0 =pz-. 

and 

(E + m)x =p</>. 

where the bispinors are given by 



and 



X= 



(ID 
(12) 



Then, assuming the non-relativistic approximation E~m, for low p, we obtain 
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from (12), and 

2 

(E-m)0=f^, (13) 

by substituting into (11). Using the same approximation, </>, here, also becomes yr. 
Conventionally, of course, the Schrodinger equation excludes the mass energy m from 
the total energy term E, and, in the presence of a potential energy V, (13) becomes 
modified to 

2 

Significant, here, is the origin of the factor 2 in the process which square roots the 

2 2 

expression E - m . The origin of the same factor in the derivation of spin from the 
Dirac equation, is seen in the behaviour of the anticommuting terms which result from 
this process. In fact, the two factors have precisely the same origin, for we may show, 

2 1 

following earlier treatments, that the Schrodinger equation can be used to derive the 
anomalous magnetic moment of the electron in the presence of a magnetic field B. 
Spin is purely a property of the multivariate nature of the p term, and has nothing to 
do with whether the equation used is relativistic or not. In our operator notation, the 
Schrodinger equation, whether field-free or in the presence of a field with vector 
potential A, can be written in the form, 

2mEy/= p 2 y/ 

Using a multivariate, p = -z'V + eA, we derive: 

2mEy/= (-z'V + eA) (-z'V + eA) yr 
= (-z'V + eA) (-z'V y/+ eA y/) 

= -V 2 y/- z'e (V. y/A + z'V y/x A + AN y/+ z'A x V y/) + e 2 A 2 y/ 

= -V 2 yr- z'e (V. y/A + 2A.V y/+ i y/V x A) + e 2 A 2 yr 

= -V 2 yr- z'e ( yN.A + 2A.V y/) + e 2 A 2 y/+ eB yr 

= (-z'V + eA).(-z'V + eA) y/+ eB yr 

= (-z'V + eA).(-z'V + eA) y/+ 2m (I.B 

This is the conventional form of the Schrodinger equation in a magnetic field for 
spin up. The wavefunction can be either scalar or nilpotent. Reversing the (relative) 
sign of eA for spin down, we obtain 
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2mEy/= (-z'V - eA) (-z'V - eA) yf 



= (-z'V - eA) (-z'V - eA) yf- 2m (i.B . 



It is significant that the standard derivation of the Schrodinger equation begins 
with the classical expression for kinetic energy p 12m, and that the factor 2 in this 
equation ultimately carries over into the same factor in the spin term for the electron. 
It is precisely because the Schrodinger equation is derived via a kinetic energy term 
that this factor enters into the expression for the spin, and this process is essentially 
the same as the process which, through the anticommuting quantities of the Dirac 
equation, makes (L + a/ 2) a constant of the motion. (Anticommuting operators also 
introduce the factor 2 in the Heisenberg uncertainty relation for the same reason, the 
Heisenberg term relating directly to the zero-point energy derived from the kinetic 
energy of the harmonic oscillator.) In fact, the origin of the factor 2, in all significant 
cases - classical, quantum, relativistic - is in the virial relation between kinetic and 
potential energies. 22 In principle, the kinetic energy relation is used when we consider 
a particle as an object in itself, described by a rest mass mo, undergoing a continuous 
change. The potential energy relation is used when we consider a particle within its 
'environment', with 'relativistic mass', in an equilibrium state requiring a discrete 
transition for any change. We can consider the kinetic energy relation to be concerned 
with the action side of Newton's third law, while the potential energy relation 
concerns both action and reaction. Because of the necessary relation between them, 
each of these approaches is a proper and complete expression of the conservation of 
energy. This fundamental relation, as we will show in sections 16 and 17, leads to the 
significant fact that the nilpotent wavefunctions, in principle, produce a kind of 
supersymmetry, with the supersymmetric partners not being so much realisable 
particles, as the couplings of the fermions and bosons to vacuum states. 



14 The hydrogen atom 



The nilpotent version of the Dirac equation has no problem in reproducing such 
results of standard Dirac theory as, for example, the hyperfine splitting of the energy 
levels in the hydrogen atom, which can be derived in straightforward terms by the 



usual textbook methods. 23 Assuming a potential energy of the form V(r), we can 
write: 

(Jeff + ijm)\]/= (kV(r) + zio.p + ijm)yf . 



Representing 



by yf and 



¥4, 



by y/' this equation decouples into two parts: 



(kH + ijm) y/ = (kV + ijm) y/ + z'io.p y/' 
(kH' + ijm) y/' = z'io.p y/ + (kV - ijm) y/' 



(14) 
(15) 



Collecting terms in (15) and dividing by 2ijm, we obtain 



(kH' -kV+ 2ijm) y/' = z'io.p y/ 
z'io.p 



( kH- kV\ 

1 +—^T. 

^ 2z/m J 



^ 2ijm ^ ' 
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The second term in the left-hand bracket is small, so. assuming the ^s to be 
eigenfunctions of the operator W with energy eigenvalues of E n , 



T lijm 



1- 



kE n -kV\. 



lijm J 



zia.p y/ 



Substituting in (14), we obtain 



kHy/ = 



kV + 



»o.p 

lijm 



1- 



kEn-kV\.. 



lijm J 



zzo.p 



J 



(16) 



Since 



d d d 



and 



0'iG.pX z'io.p) = p 

2 i dV 

(z'io.p) V(r)( z'io.p) = Vp --^(o.rX o.p) 



= V P -- ¥ (r.p)-- ¥ (o.rxp) 



Equation (16) now becomes 



kHy/ = 



kV + 



lijm 



kE n -kV } 
lijm 



i kdV 



1 A:3V 



4m 



^¥ r - p+ 4^7¥ (o - rxp) r - (17) 



Using S = a/2 and L = r x p, and taking £„ - V as the nonrelativistic kinetic energy 
term p 2 /lm, equation (17) becomes 



kfTy/ = 



kV + 



lijm 



1- 



^pT| k BVd_ 1 kdV 



1m J j Am 2 dr dr + lm z r dr 



S.L 



(18) 



The fourth term in the bracket represents k times the interaction energy between 
spin and orbital moments. With a Coulomb potential due to an isolated nucleus of 
charge Ze, 



V(r) = 



Ze L 
4ft£ r 



this energy becomes 



1 kdV „ Ze 2 1 „ 
S.L = — 2 3 S.L 



2m r dr 



In the standard theory, using the hydrogenic wavefunctions as derived from the 
nonrelativistic Schrodinger theory, the mean value of 1/r 3 is calculated as 
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ryi 3 6 

f Hyr) ¥ MV= ( ^ £o)Vl(l + m+l) 

oo 

where m r is the reduced mass (~ m), and / > 1, because the term vanishes when / = 
because L.S vanishes. The average interaction energy between L and S then becomes 

ry 2 ryZ 3 6 ^4 8 

ze L m r e z me 

8^ m 2 (4^)V/(/ + Vi)(l + 1) = 2{%ne ) 4 nl{l + + 1) S ' ' 

With the Rydberg constant for infinite mass given by, 

4 

me 

and the fine structure constant, 



we obtain 

Z^mg 8 _ InR^dz 4 

2(Sx£ ) 4 n 3 l(l + Vi){l + 1) SL ~ n 3 /(/ + 1 /2)(/+l) SL ' 

the energy equation for the hyperfine splitting. 

15 The origin of the factor 2 in spin states 

We have already suggested that all the important factors of 2 (or Vi) in classical 
physics, relativity and quantum physics, can always be attributed to whether one is 
using potential or kinetic energies, or using just the action side, or both action and 
reaction sides, of Newton's third law, and that this is a product of whether one is 
using continuous or discrete solutions, or changing or fixed ones. Through a series of 
parallelisms and derivations, we can show that this factor 2 comes from the symmetry 
between the action of an object and the reaction of its environment (whether material 
or vacuum), and that, while the object taken on its own would show changing 
behaviour, leading to an integration resulting in a factor Vi (kinetic energy), a 
conservative 'system' of object plus environment would show unchanging behaviour, 
expressible in terms of a potential energy term, which is twice the kinetic energy. 

The deeper aspect to this argument is that it makes sense of the boson / fermion 
distinction in a fundamental way, as well as supersymmetry, vacuum polarization, 
pair production, renormalization, spin Vi, and so on, because the halving of energy in 
'isolating' the fermion from its vacuum 'environment' is the same process as 
mathematically square-rooting the quantum operator via the Dirac equation. The 
work of Bell et al. 17 on producing integral spins from half-integral spin electrons 
using the Berry phase is a particular case of this, and it is through this kind of case 
that we may extend the principle in the direction of supersymmetry. 

The act of defining a rest mass also defines an isolated object, and one cannot 
define kinetic energy in terms of anything but this rest mass. If, however, we take a 
relativistic mass, we are already incorporating the effects of the environment. The 
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most obvious instance is that of the photon. The photon has no rest mass, only a 
relativistic mass; mc 2 for a photon behaves exactly like a classical potential energy 
term, as well as having the exact form of a potential energy for a body of mass m and 
speed c. A particular instance is the application of a material gas analogy to an photon 
gas in producing radiation pressure pc 13. Action and reaction occurs in this instance 
because the doubling of the value of the energy term comes from the doubling of the 
momentum produced by the rebound of the molecules / photons from the walls of the 
container - a classic two-step process, like the two-way speed of light. 

The energy involved in both material and photon gas pressure derivations is 
clearly a potential energy term (the material gas energy has to be halved to relate the 
kinetic energy of the molecules to temperature), and its double nature is derived from 
the two-way process which it involves, which is the same thing as saying that it is 
Newton's action and reaction. The same thing happens with radiation reaction, which 
produces a 'mysterious' doubling of energy hv I 2 to hv in many cases, and also 
zitterbewegung for the electron, which is interpreted as a switching between two 
states. Feynman and Wheeler also produce a doubling of the contribution of the 
retarded wave in electromagnetic theory, at the expense of the advanced wave, by 
assuming that the vacuum behaves as a perfect absorber and reradiator of radiation. In 
principle, this seems to be equivalent to assuming a filled vacuum for advanced 
waves (equivalent to Dirac's filled vacuum for antimatter), and relates to previously 
stated ideas that continuity of mass-energy in the vacuum is related to the 
unidirectionality of time. 4 ' 5 ' 6 

In the case of material bodies, an old proof of Newton's of the mv 2 lr law for 
centripetal force, and hence of the formula mv for orbital potential energy, had the 
satellite object being 'reflected' off the circle of the orbit, first in a square formation, 
and then in a polygon with an increasing number of sides, becoming, in the limiting 
case, a circle. Here, the momentum-doubling action and reaction produce the 
potential energy formula, and demonstrate the relation between the conservation laws 
of linear and angular momentum. 

Interestingly, though light in free space has velocity c, and, through special 
relativity, has no rest mass, and therefore no kinetic energy, as soon as you apply a 
gravitational field, the light 'slows down', and, at least behaves as though it can be 
treated as a particle with kinetic energy in the field. Immediately, then, you can apply 
a kinetic energy-type equation, and derive the double gravitational bending 
(ultimately from the 2 in mc 2 /2) using the very simple mc 2 /2 = GMm (e - l)/R, where 
e is the eccentricity of the hyperbolic orbit. 

The result of all these cases is that kinetic energy variation may be thought of as 
continuous, but starting from a discrete state; potential energy variation, on the other 
hand, is a discrete variation, starting from a continuous state. Each creates the 
opposite in its variation from itself. Kinetic energy and potential energy create each 
other, in the same way as they are related by a numerical relationship. 

The existence of two conservation of energy approaches thus has very profound 
implications, and appears to arise from a very deep stratum in physics. Kinetic energy 
may be associated with rest mass, because it cannot be defined without it - one could 
consider light 'slowing down' in a gravitational field as effectively equivalent to 
adopting a rest mass, and, of course, photons do acquire 'effective masses' in 
condensed matter. Potential energy is associated with 'relativistic' mass because the 
latter is defined through a potential energy-type term (E = mc ), light in free space 
being the extreme case, with no kinetic energy / rest mass, and 100 per cent potential 
energy / relativistic mass. The description, in addition, seems to fit in with the halving 
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that goes on, for a material particle, when we expand its relativistic mass-energy term 
(mc 2 ) to find its kine 
conservation equation 



2 2 

(mc ) to find its kinetic energy (mv 12). Thus, if we take the relativistic energy 



E 2 - p 2 - m 2 = , 

we can take regard this as a 'relativistic' mass (potential energy equation of the form 
E = mc 2 (treating at one go the particle interacting with its environment), and proceed 
to quantise to a Klein-Gordon equation, with integral spin. Alternatively, we can 
separate out the kinetic energy term using the rest mass m . From 



E = mo c 



v c 



We take the square root, and obtain 

2 

2 mv 

E = moc + ~7T~ + 



In this case, we are looking at the 'isolated' system of the electron itself. The 
Schrodinger equation, of course, arises from this approach, quantising mv 2 l2 in the 
form p 2 l2m. Now, as we have shown, using a multivariate form of the momentum 
operator, p = -z'V + eA, the Schrodinger equation produces the magnetic moment of 
the electron, with the required half-integral value of spin, the Vi coming from the term 

2 2 

mv 12 or p 12m; it is also effectively a limiting approximation to the bispinor form of 
the Dirac equation. In principle, this means that the spin Vi term that arises from the 
Dirac equation has nothing to do with the fact that the equation is relativistic, but 
arises from the fundamentally multivariate nature of its use of the momentum 
operator (equivalent to the use of Pauli matrices). In fact, it comes from the very act 
of square rooting the energy equation in the same way as that operation produces 
mv 12 in the relativistic expansion. The Vi is, in essence, a statement of the act of 
square-rooting, which is exactly what happens when we split into two nilpotents; 
the Vi in the Schrodinger approximation is a manifestation of this which we can trace 
through the Vi in the relativistic binomial approximation. 



16 The fermion and its 'environment' 



Our Dirac equation 



( 3 ^ 
iter. + i V+ ijm yf = , 

V dt J 



was, of course, obtained by direct square-rooting of the equation 



E 2 -p 2 -m 2 = 



and leads to a column vector of four terms: 



(kE+iip + ijm) e - i(E '- p - r) 
(kE-iiv + ijm) e - i(Et - p - r) 
(-kE + ii p + ij m) e - i(Et - p - r) 
(-kE-iiv + ijm)e i(Et - p - r) , 
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representing a single quantum state. We have shown that a spin 1 boson creation 
operator is of the form 

(kE + ii p + ij m) (-kE + ii p + ij m) 
(kE - ii p + ij m) (-kE - ii p + ij m) 
(-kE + ii p + ij m) (kE + ii p + ij m) 
(-kE - ii p + ij m) (kE - ii p + ij m) 

while a spin boson can be represented by 

(kE + ii p + ij m) (-kE - ii p + ij m) 
(kE - ii p + ij m) (-kE + ii p + ij m) 
(-kE + ii p + ij m) (kE - ii p + ij m) 
(-kE - ii p + ij m) (kE + ii p + ij m) , 

each multiplied by the usual exponential form in creating the wavefunction. While the 
fermion wavefunction is effectively a nilpotent or square root of 0, the boson 
wavefunction is a product of two nilpotents, each not nilpotent to the other. The 
multiplications here are scalar multiplications of a 4-component bra vector (composed 
of the left-hand brackets), representing the particle states, and a ket vector (composed 
of the right-hand brackets), representing the antiparticle states. 

Though we generally consider fermions to follow the Dirac equation, derived, as 
we have seen, from kinetic energy, and bosons to follow the Klein-Gordon equation, 
derived ultimately from potential energy, when we look at real situations something 
remarkable happens. Fermions with spin Vi become spin 1 particles when taken in 
conjunction with their 'environment' - the Jahn-Teller effect (via the Berry phase) 
and Aharanov-Bohm effect are examples. A considerable number of phenomena are 
already associated with the Berry phase, and there are no doubt a large number of 
others waiting to be discovered. 

Bosons with spin or 1 become spin Vi or 3/2 by the reverse process - this is 
what happens when the gluons interact with quark-gluon plasma to deliver spin Vi or 
3/2 to a baryon. It is almost certainly a universal principle that fermions / bosons 
always produce a 'reaction' within their environment, which couples them to the 
appropriate wavefunction-changing term, so that the potential / kinetic energy relation 
can be maintained at the same time as its opposite. We can relate this to the whole 
process of renormalization, which produces an infinite chain of such couplings 
through the vacuum. (The related Feynman- Wheeler process of vacuum absorption of 
radiation, of course, also reduces the infinite electron self-energy to a finite mass.) It 
is also of significance here that the vacuum wavefunctions we have previously 
defined are of the complementary forms: (-kE + ii p + ij m) for fermions, and (kE + ii 
p + ij m) for antifermions. 

17 Supersymmetry and renormalisation 

We have seen that Vi spin is always due to kinetic energy (from continuous variation) 
and we can consider it as the isolated fermion; the unit spin comes from potential 
energy (stable state) and represents either a boson with two nilpotents (which are not 
nilpotent to each other) or the equivalent, which is a fermion interacting with its 
environment - and manifesting Berry phase, radiation reaction, relativistic correction, 



31 



Thomas precession, or whatever else is needed to produce a 'conjugate' spin state in 
its 'environment' (vacuum). It is a question of whether we are treating the action half 
of Newton's third law, or the action and reaction pair. The 'supersymmetric' partner 
seemingly comes from this choice of fermion or fermion plus environment approach. 

The nilpotent operators defined for fermions are also supersymmetry operators, 
which produce the supersymmetric partner in the particle itself. The advantage of the 
supersymmetry concept is that boson contributions and fermion contributions are of 
opposite sign (with the operators having opposite signs of E) and automatically 
cancel. The particular advantage of this realisation is that there are no extra 
supersymmetric partners yet to be discovered. 

The Jahn-Teller effect throws up the Berry Phase phenomenon because it is 
treated semi-classically by making an adiabatic approximation (AA). Assuming the 
time scale for nuclear motions is very much greater than the time scale for electronic 
transitions, we get a total wave function made up of three factors: 

1) A phase factor 

2) A factor associated with the electronic coordinates (q) 

3) A factor associated with the nuclear coordinates (Q) 

The Jahn-Teller effect couples the {Q} and {q} motions so that different parts of the 
total wave function change sign in a coordinated manner to preserve the single- 
valuedness of the total wave function. Neither the nuclear nor the electronic 
wavefunction are single-valued by themselves, though the total wavefunction is. 

Thus, a sign change occurs when a fermion rotates through 2n, but if the 
relationship between fermion and 'the rest of the universe' is like the total wave 
function of the Jahn-Teller effect then there is a preservation of single-valuedness. 
The Jahn-Teller effect is also the spontaneous breaking of a point-group symmetry, 
the electronic part of the Hamiltonian producing a Mexican hat potential upon which 
the nuclei move. This adiabatic potential is the basis of the symmetry breaking. It is 
the Jahn-Teller effect that couples the nuclear coordinates to the electronic 
coordinates, and it may be that the same thing be going on in continuous symmetries 
to couple electrons to the vacuum. (Certain vibronic - vibrational / electronic - 
interactions take up continuous group symmetries to first order in the Jahn-Teller 
effect, and this may have parallels with the Higgs mechanism.) 

If this is true, the coupling of the vacuum to electrons could be generating 
'boson-images' and vice versa. Presumably the required loop diagrams that lead to 
renormalisation (cancellation of fermion with boson loops) would then have an 
explanation without us needing to postulate the existence of extra boson/fermion 
equivalents (photon/photino, etc.). 'Supersymmetry' may be part of a much more 
general pattern. Bosons and fermions seem to require 'partner states' as much as 
potential and kinetic energy are needed to fully describe conservation. 

The Q generator for supersymmetry is simply the term (kE + zip + ijm). To 
convert bosons to fermions, or antifermions to bosons, we multiply by (kE + zip + 
ijm) (the p could, of course, be + or -). To convert bosons to antifermions, or 
fermions to bosons, we multiply by (-kE + z'ip + ijm). These are Q and Q\, and 
(taking the normalisation that we have previously introduced for the vacuum 
operator) the anticommutator becomes effectively E, the Hamiltonian. Of course, if 
we write out the operators in full, they will be respectively four-term bra and ket 
vectors, with the E and p values going through the usual cycle of + and - values. 

Using the 'environment' principle, we can imagine an infinite series of 
interacting terms of the form (kE + z'ip + ijm), (kE + zip + ijm) (-kE + zip + ijm), (kE 
+ z'&p + ijm) (-kE + z'ip + ijm)( kE + z'ip + ijm), (kE + z'ip + ijm) (-kE + z'ip + ijm)( kE 
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+ zip + ijm) (-kE + zip + ijm), etc., where selection of the appropriate terms will lead 
to a cancellation of the boson and fermion loops of opposite sign. The (kE + zip + 
ijm) and (-kE + zip + ijm) or bra and ket factors are an expression of the behaviour of 
the vacuum. To convert bosons to fermions we multiply by (kE + zip + ijm) (the p 
could, of course, be + or -). To convert bosons to antifermions, or fermions to 
bosons, we multiply by (-kE + z'ip + ijm). 

In effect, the vacuum state must act like a 'mirror image' to the fermion. 

a | 0> = (kE + z'ip + ijm) k (kE + z'ip + ijm) 

not only looks like a 'mirror image' when you write it down, but is the expression of 
part of an infinite regression of images of the form 

(kE + z'ip + ijm) k (kE + z'ip + ijm) k (kE + z'ip + ijm) k (kE + z'ip + ijm) ... 

Though it is possible to construct a vacuum state that these operators act upon, one 
gets a vacuum state that is dependent on the operator that is acting on it. Thus, the 
vacuum state of the state (kE\ + z'ipi + ijm) turns out to be k (kE\ + z'ipi + ijm). In 
addition, however, 

(kE + z'ip + ijm) k (kE + z'ip + ijm) k (kE + z'ip + ijm) k (kE + z'ip + ijm) ... 

is the same as 

(kE + z'ip + ijm) (-kE + z'ip + ijm) (kE + z'ip + ijm)( -kE + z'ip + ijm) ... 

The infinite series of creation acts on vacuum is the same as the mechanism for 
fermion / boson / fermion / boson ... creation (supers ymmetry, renormalization), and 
is only true if the series is infinite, because each bracket has to be postmultiplied by k 
to make it work. The spin terms p must, of course, be of the same sign to produce spin 
1 bosons - spin 0, as in the case of the mass -generating Higgs, would break the 
sequence. The 'reflection' process also implies an infinite range of virtual E values in 
vacuum (though, presumably, only of one sign in the ground state). 

Here, then, we may even have an explanation of the mechanism of 
renormalisation. It is also worth noting that the Klein-Gordon equation automatically 
applies to fermions (because it simply involves pre-multiplication of zero by a 
nilpotent differential operator), while the Dirac equation applies to 'spin 1' particles 
(in the particle + environment sense - like a method of images), leading to the 
Aharonov-Bohm effect, the Jahn-Teller effect, zitterbewegung, etc. With the 
'environment' principle we don't need to suppose we have an extra set of bosons or 
fermions - the coupling to the vacuum is automatic in a series of entangled states. 

In splitting away a fermion from the 'system' we introduce a coupling as a 
mathematical description of being able to split it away; the coupling is to the rest of 
the universe to make the total wavefunction single-valued. Fermions by themselves 
cannot have single-valued wavefunctions, if this total wavefunction is to be single- 
valued. It may be possible, as we have supposed, to predict effects of the Jahn-Teller, 
Aharonov-Bohm type, where the fermion acquires a single-valued wavefunction by 
interaction with its 'environment'. There must be other effects of this type which can 
be predicted and detected - maybe many others. 
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The reverse effect must also be possible, of bosons coupling to an 'environment' 
to produce fermion-like states. The most immediately likely possibility is the 
behaviour of gluons in producing the baryon total spin of Vi or 3/2. The six- 
component baryon wavefunction has something equivalent to (kE ± iip x + ijm) (kE ± 
iipy + ijm) {kE ± iip z + ijm), and the p x , p y , p z and ± represent the six degrees of 
freedom for p. We can imagine the p rotating through the three spatial positions 
leaving terms like (kE + Up + ijm) (kE + ijm) (kE + ijm); (kE + ijm) (kE ± Up + ijm) 
(kE + ijm); etc., which effectively produce something like a fermion combined with 
Bose-Einstein condensate. This is correct, but how do the gluons do it? Presumably, 
in effect, by the gluon 'transferring' the p between one (kE + ijm) and another - to 
use crudely mechanistic terms - becoming a boson of spin 1 with an effective 
contribution from the 'environment' due to the gluon sea making it transfer spin Vi. 

Another aspect of the process is that dimensionality, in general, introduces two 
orders of meaning in a parameter - of the value (as in length/time or charge/mass), 
and of the squared value (as in Pythagorean/vector addition of space dimensions, or 
space and time, or energy and momentum, or charges/masses 'interacting' to produce 
forces). In a sense we are doing this with fermion and boson wavefunctions, one type 
being a 'square root' of the other. 

18 The creation of the Dirac state 

We have now reached the point where we can suggest how the algebra is 
involved in the creation of the Dirac state. We have four fundamental parameters, 
time, space, mass, charge, which require the respective algebraic units i, i, j, k, 1, i,j, 
k. However, the fundamental units of the algebra which they entail is created from an 
anticommuting pentad, with specific rules concerning the two 3-dimensional 
arrangements, i, j, k, i,j, k? To reduce the 8 pure units to 5 composite ones, we need 
to map one of them onto the other three parameters. So, we map charge onto space, 
time and mass, w, s, e, each taking one of the parameters (though we could no doubt 
also obtain something useful by mapping space onto time, charge and mass). This 
gives us an algebra for the composites of the form ki, ii, zj, ii,j, which is, of course, 
the Dirac algebra (though we usually there multiply by an extra i purely for 
convenience). 

But what have we done physically? We have, by putting charge components onto 
time, space and mass, actually introduced quantization (for charge is, of course, itself 
quantized), and since the charges are conserved quantities, we have also created a 
quantum state with fixed E-p-m. The system, of course, is already relativistic because 
our original units were, and the exponential term serves to preserve the required 
variation in time and space. £-p-m is the direct result of mapping conserved unit 
charges onto time, space and mass. The concept of rest mass is only created by the act 
of 'quantization'. 

Charge, being irrotational, because of the separate conservations of the weak, 
strong and electric charges, should not rotate the values of w-s-e, but the actual 
mapping is arbitrary. Hence we have to have rules for 'charge accommodation'. The 
k, i,j in the Dirac equation are, in principle, the charges (and this, in a sense, justifies 
the Dirac equation's otherwise arbitrary interpretation of -E states as antiparticles), 
and the £-p-m Dirac state is actually created by w-s-e. This is a direct link between 
E-\)-m and w-s-e, in addition to justifying the fundamental nature of the time, space, 
mass, charge paradigm, and explaining in physical terms why we have such a thing as 
a Dirac state. 
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E-p-m and w-s-e reflect different ways of combining the 8 parts to make 5 (we 
may refer to this process as 'compactification' or 'folding'), and so one might expect 
different mathematical structures as we make the combination E-p-m-like or w-s-e- 
like. Doing the 'compactification' in two different ways doesn't mean that the 
resulting composite notions will map directly onto each other. E and p are fixed in a 
Dirac state, but they have an infinite range of possible values. This is not true of w 
and s, and we cannot necessarily infer that the value of £7p will reflect that of wis (or 
even its presence). £-p-m is certainly easier to generate than w-s-e, but one can see 
that the process of 'folding' the charge quaternions onto time, space, mass can be said 
to make one charge time-like, another space-like, and another mass-like. It may be, 
however, that the 'folding' process may be of the charge quaternions onto the Dirac 
E-\i-m, rather than the original time, space, mass. 

£-p-m is discrete, being made so in the process of 'compactifying' the charge. If 
imaginary numbers create a discontinuous U(5), 24 then this is because in E-p-m we 
have a real quantized mass and an imaginary quantized E. Imaginary i is not only 
necessary to complete the Dirac pentad, it is also necessary to fermions because E I w 
is necessary to fermions. It is also interesting that the group is U(5) because the 
gravity operator ('mass') has been automatically included in the process of 
'compactification'. Imaginary / must introduce discreteness in the square-rooting 
process which produces iE, and this links up with the fermion / boson distinction - 
square-rooting the energy equation introduces discreteness (the fermion, with rest 
mass) in a continuous variation (the kinetic energy). 
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